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Abstract 

The ground state of an electrorheological (ER) fluid has been studied based 
on our recently proposed dipole-induced dipole (DID) model. We obtained an 
analytic expression of the interaction between chains of particles which are of 
the same or different dielectric constants. The effects of dielectric constants on 
the structure formation in monodisperse and polydisperse electrorheological 
fluids are studied in a wide range of dielectric contrasts between the particles 
and the base fluid. Our results showed that the established body-centered 
tetragonal ground state in monodisperse ER fluids may become unstable due 
to a polydispersity in the particle dielectric constants. While our results agree 
with that of the fully multipole theory, the DID model is much simpler, which 
offers a basis for computer simulations in polydisperse ER fluids. 

PACS Numbers: 83.80.Gv, 82.70.Dd, 41.20.-q 
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I. INTRODUCTION 



The study of structure formation of electrorheological (ER) fluids has attracted increasing 
interest in recent years for its fundamental and technological values. Upon the application 
of an external electric field, the suspended particles in an ER fluid aggregate into chains 
and then columns parallel to the field, and drastically change the rheology of the suspension 
[1]. Tao and his coworker first suggested the existence of microcrystalline structures inside 
the columns and identified its ground state to be a body-centered tetragonal (bet) lattice 
[2]. These authors developed an analytic method based on a point-dipole (PD) approxi- 
mation, i.e., treating the dielectric spheres as point dipoles interacting with one another. 
The possible ground state is the configuration that minimizes the dipole interaction energy 
and consequently the total Coulomb energy. The idea was soon confirmed by computer 
simulations and experiments [3,4]. The discovery of this property is not only helpful to 
understand further the ER response to the external field, but also offers a new technique to 
form mesocrystals with unique photonic properties [5,6]. 

After many efforts to reveal the details of structure formation in an ER fluid, it has 
been known that what structure will be formed in such a system actually depends on many 
factors, such as volume fractions [7], the size distribution of particles [8,9], the field frequency 
[10], etc. Combining the electrorheological and magnetorheological effects [5,11] or applying 
a rotating electric field [12] was also found to cause a structure transition from the bet to 
face-centered cubic (fee) lattices [13]. 

The influence of the dielectric constant is another interesting topic since dielectric mis- 
match between the particles and the fluid is widely accepted as the main reason for the 
ER phenomenon. Davis found that the bet, fee and hexagonal closed packed (hep) lattices 
degenerate when the particle permittivity (e p ) is much larger than that of the base fluid {ef) 
[14]. Clercx and Bossis [15] developed a fully multipolar treatment and compared the bet 
and fee lattices for some values of e p /e/. Their results are in agreement with Davis. Lukkari- 
nen and Kaski studied both the free energies of many types of lattices and dynamical effects 
in an ER fluid containing particles with the dielectric constant being greater and/or smaller 
than that of the fluid [16,17]. However, all of the previous work has either limited to the 
extreme case (e p ^> ef) or been too complicated to be adopted in further studies. And the 
understanding of the way dielectric constant affects structure formation is far from being 
sufficient. 

Apart from theoretical work, computer simulation is another effective method to gain 
insight into ER effects. The point-dipole approximation is routinely adopted in simulation 
[18] because of its simplicity. Since many-body and multipolar interactions are neglected 
in this approximation, the predicted strength of ER effects is of an order lower than the 
experimental results . On the other hand, the accurate theoretical models are usually too 
complex to use in dynamic simulation of ER fluids. Hence, a model which is easy to use but 
beyond the point-dipole results is needed in computer simulations. 

In this paper, we use a dipole-induced-dipole (DID) model [19] proposed by one of the 
authors to study the ground state of ER fluids in a wide range of the permittivity ratios: 
< e p /ef < oo. This model accounts for the multipolar interaction but is significantly 
simpler compared with existing multipolar theories [15,20,21], therefore it can serve as a 
candidate for computer simulation instead of the traditional point-dipole approximation. 
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Some computation has been carried out such as the calculation of the interacting force 
between particles of different sizes and various dielectric constants and the simulation of the 
athermal aggregation of particles in ER fluids, both in uniaxial and rotating fields [19,22,23]. 
The purpose of this paper is to use this model to deal with structure formation in ER fluids 
and to offer some theoretical prediction as instructions to further dynamic simulations of ER 
fluids. It is found the ground state of ER fluids may vary with the dielectric constant, and 
the critical dielectric contrast is estimated. Furthermore, we obtain the interaction between 
two chains containing particles of various dielectric constants and apply it to polydisperse 
ER fluids where the particles have the same size but different permittivities. The effects of 
polydispersity on structure formation are investigated. 



II. INTERACTION BETWEEN TWO PARTICLES 

We start with briefly reviewing the DID model and then applying it to deal with the 
interaction energy between two particles. The DID model is defined from a multiple image 
method [19] which is based on a generalization of the image method to dielectric spheres. 
First consider a simple situation in which a point dipole p is placed at a distance r from 
the center of a sphere. The orientation of the dipole is perpendicular to the line joining the 
dipole and the center of the sphere. If the sphere is conducting, the induced image dipole is 
exactly given by p' = —p(a/r) 3 and at a distance r' = a 2 /r from the center. Generalizing this 
result to a sphere of dielectric constant e p placed in a host medium e/ reads p' = — j3p(a/r) 3 
with 3 as the dipolar factor 3 = £p ~ £/ . If the orientation of the point dipole is parallel to 

€p-\-Z€f 

the axis, then p' = 28p(a/r) 3 . In the limit f3 — > 1, the above results reduce to the conducting 
sphere case. 

Then consider a pair of dielectric spheres A and B, of the same radius a but different 
dielectric constants t v \ and e P 2, separated by a distance r, in a base fluid of a dielectric 
constant e/. Upon the application of an electric field E , the induced dipole moments in the 
individual spheres are respectively given by 

Pao = e f pE a 3 ,p b0 = tfp'Eoa 3 (1) 

with the dipole factors 3 = £pl ~ <:/ and 3' = £p2 ~ £/ . The initial dipole moment p a0 induces 
an image dipole p b i in sphere B, while pbi induces another image dipole p a 2 in sphere A. As 
a result an infinite series of dipoles p a o,p a i,Pa2, ■■■ are formed inside sphere A and the total 
dipole moments for transverse and longitudinal fields are respectively given by 

p aL = sinh 3 a f [ P'°( 2 fl W ~W) w - 1 + Vma\WnWT-\ (2) 
~J (sinhna + sinh(n — l)a) 3 (rsinhna) 3 ' 

p aT = sinh 3 a f \ P^-W't-f 3 '^ 1 + P^ 3 (-Pn-PT~\ ^ (3) 
^ (sinh no; + sinh(n — l)o;) 3 (rsinhna) 3 

The parameter a in Eq. (1) and (2) satisfies 

r 2 _ 2a 2 

cosh a = — — — . (4) 
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Similar expressions for the total dipole moments inside sphere B, namely pbL and PbT, can 
be obtained by interchanging (3 and /3'. 

We should remark that the present generalization is only approximate because there is 
no image method for a dielectric sphere. An integral equation approach was proposed to 
examine the validity of the multiple image method in Ref. [19]. It was shown that Eq. (2) 
and (3) can produce good results at high dielectric contrast j3 — > 1 as expected, and the 
model keeps reasonable even in the low contrast case [19,22]. 

How many terms are retained in the analytic multiple image result determines different 
models. With only the n — 1 term in the series of Eq. (2) and (3), it reduces to the 
traditional point-dipole approximation, while the DID model is defined by retaining the 
first two terms (n — 1 and n — 2). Taking part of the multipole effects into account, the 
DID model is generally better than the PD results, especially in polydisperse systems, and 
can be chosen as an effective approximation to calculate multipolar interaction [19,22,23]. 

The incorporation of multipole effects leads to one of the obvious differences between the 
PD and DID model: the symmetry of negative and positive (3s will be broken in the latter, 
which is illustrated in Fig. 1. This asymmetry of the dipolar factor will also be exhibited 
in the interaction energy between particles, as shown below, and finally affect the structure 
formation in ER fluids. It implies that different structures may form when the dielectric 
constant of particles is smaller than that of the host fluid. 

Now we begin to apply the DID model to deal with the interaction energy between two 
particles. The electrostatic energy of the two particles upon the application of an electric 
field E at an arbitrary angle 9 (Fig. 2a) is given by [24] 

E 

w = —£[{PaL + Pb L ) cos 2 9 + ( PaT + p bT ) sin 2 9}. (5) 

Substituting the dipole moment expressions into Eq. (5) and letting A = f3'//3, po = p a o, we 
can get the energy using the DID approximation as 



1 + A ,_, . , pi 3 cos 2 9 — 1. A + A 2 r pi 

W = —E oPo + A(-— + — — - [--^PfJ. 3 cos 2 9 + 1 

€f r 6 2 e/cr 



+ A 2 [-^L 2 z/(9cos 2 £-l)] 

= Wf + W + Wi + W 2 , (6) 
where the dimensionless parameters have been defined: 

, , 2 1 t) 3 > " = (^) 3 > r ' = 

r z — 1 r z — 2 a 



^ = (3^T) 3 > -=(3^) 3 > ' = (7) 



The first term Wf in Eq. (6) is the energy of the individual particles by the applied field, 
which is independent of the relative position of the particles, while the remaining terms 
correspond to the interaction energy which determines the formation of ground states [2]. 

Note that the term Wq is, if A = 1, exactly the result of Tao et a/.'s work [2] , in which 
they dealt with the monodisperse case in the use of a PD approximation and concluded that 
a bet lattice will form as the ground state. So what we are interested in is the last two terms 
which describe the DID correction to the total interparticle energy. Adding the hard-core 
repulsion which keeps the particles from coalescing (r' > 2), W\ and W 2 can be rewritten as 
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.A + A 2 ' 3cos 2 0+l \ + \ 2 ' r , 3 3,1 

Wl = -0_ £ „„_, —p^- = -ft— £ a^l4 + -P^]- (8) 

W„ = -0 2 \ 2 ± a„J^^± = -p 2 \ 2 ± a^[8 + IpAjij, (9) 

n=9 ' n=9 ' 4 ^ ' 

Here we have taken as the energy unit. J2' n =i(---) means the summation over n = 

i,i + 2,i + 4:, the coefficient a; = 3(3+1) -., (3+ '~ 1} , and p = [(x z - x b f + (y a - y b ) 2 ] 1/2 = 

\Jr 2 — (z a — Zb) 2 ((xi, Zi) with i = a,b being the coordinates of the particles and z axis is 
chosen as along the field direction). Obviously, the DID correction cannot be ignored when 
(5 is large, while the term W\ being proportional to f3 instead of \/3\ shows the asymmetry 
of (3 in the DID model. 

III. MONODISPERSE SYSTEMS 

Now we study a monodisperse ER fluid composed of particles with a dielectric constant 
e p , which can be larger (f3 > 0) or smaller (f3 < 0) than that of the fluid e/. It has been 
mentioned that the interacting particles in a ER fluid will first form chains between electrode 
plates, then the chains aggregate into columns containing microstructures. Combined with 
the images, the chains can be treated as infinite if the distance between two plates is large 
enough [2]. The interaction energy per particle is divided into two parts: one is from the 
self-energy of an infinite chain, i.e., the interaction energy between the particles belonging 
to the same chain; the other from interaction between different chains [2]. 

A. u s , the self- energy 

Consider an infinite chain containing particles of radius a and dielectric constant e p at 
ij = 2a jz (j — 0, ±1, ±2, ...) (Fig. 2b). The self-energy of the chain is given by 

« S = ^E Wifo) = \ E Wi{r = 2aj, 6 = 0). (10) 

Wi refers to the interaction energy between two particles. Substituting the DID terms in 
the particle interaction, W\ and W2, into Eq. (10), we get 

' G U-6 

< = ~P E oS=2 C(n) = -0.149449/5, (11) 

n=6 

' fl n-9 

u 2 = -P 2 E TSrCW = -0.125047/3 2 , (12) 

n=9 2 2 

where ((n) is the Zeta function defined as = Yl™=i jr- Again, the DID correction is 
comparable to the PD result [2], —0.300514, when (3 is large. The total self-energy u s in a 
DID model is the sum of u\ , uf, and that from PD assumption 

u s = -0.300514 - 0.149449/? - 0.125047/3 2 , (13) 

which represents a correction to the established PD results. 



5 



B. u l , the interaction energy 

The interchain energy between two parallel infinite chains, with vertical shift z and 
separated by distance p, is given by ju l (p, z), where u l (p,z) is the interaction between one 
dielectric particle at r p = p + zz near the infinite chain in which particles locate at ij = 2ajz 
(j=0,±l,±2,...) (Fig. 2c), 

= EWifo-r,). (14) 

3 

Hence we have the DID correction of u l 

' 3 d +oc 1 

u\(p, z) — —3 V an- 6 [4 + —p— I V ^-r , ,. 01 ;o (15) 

and 

ui(p,z) = -3 2 Y an_ 9 [8 + -p^-] Y — r^r. (16) 

2 ^ ^ n=9 ~ n V ^oc [P 2 + (*- 2«j) 2 ] n/2 1 ^ 

Following a Fourier series technique proposed by Tao et al. [2] we expand u\ and i4 into 



with 



51 = y~^(— )~ r " [(3 + 5n)if n-i (pa;) — Spw-fTn+i (pa;) — ?>puiK n -3 (pa;)] cos( ), (19) 

s=1 2 2 2 2 ' v a 

5 2 = £(— ^[fa + 7n)K^i(pu) - VpuKn+iipuj) - 9pujKn^(pcu)] cos(— ), (20) 

~[ 2 2 2 2 a 

a; = sn/a. (21) 

T(x) and Ki(x) in the above equations are the gamma function and i-th order modified 
Bessel function, respectively. The sums above can be easily evaluated numerically, and the 
expression of u' l (p, z) is given by u l (p, z) = u l Q (p, z) + u\(p, z) +u l 2 (p, z), where u l (p, z) is the 
PD result as [2] 

u »(p, Z ) = \ £ 2 **s 2 K (^) «*(—). (22) 

8 = 1 a a 
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C. Possible ground states 



The total interaction energy per particle in a certain configuration is [11] 

1 ' 

u = u s + -Y,u l (Pk,Zk), (23) 
1 k 

where J2'k denotes the summation over all chains labeled k except the one containing the 
considered particle. Since the self-energy of a chain is independent of the structure, what 
really affects the energy difference between various lattices is u l (p,z). Fig. 3 shows the 
dependence of u l on the shift z for different values of p when f3 > 0. The results from 
the PD approximation are also plotted for comparison. The interaction may be either 
attractive or repulsive depending on the shift z, and the range of z in which two chains 
attract each other enlarges quickly when (3 increases. This implies a tendency to form a 
more close packed structure than the bet lattice with increasing dielectric contract. The 
most possible candidate is the fee lattice, which is nearest to the bet one in energy [2]. An 
estimation including the nearest and the next-nearest neighboring chains can give the energy 
gap between the two lattices 

Au = Uf cc — u bct = u l (p = 2a, z = 0) — 2u l (p = V6a, z — 0) 

= 0.0110 -0.0230/3 + 0.0127/3 2 , (24) 

which decreases when (3 increases from to 1. For the limit case f3 — 1 (e p >■ e/), Au ~ 6.6 x 
10~ 4 and the two phases can be regarded as degenerate, in agreement with the conclusion 
of previous work [14,15]. 

Reversed situation happens when f3 < (Fig. 4 ). Now it's the repulsive region that 
enlarges and the system tends to form a looser structure as (3 approaches the negative limit 
— \ (e p <C €f). The loosest structure is certainly the one composed of separate chains. 
Calculating the difference between the bet lattice and separate chains, we get 

Au = 2u\p = Via, z = a)+ 2u\p = VGa, z = 0) 

= -0.081 - 0.300/3 - 0.155/3 2 . (25) 

When f3 = —0.323 a transition from the bet lattice to separate chains happens and the latter 
becomes more and more stable while the absolute value of f3 increases. 

Finally we want to make a comparison with Clercx et a/.'s elegant work based on a 
multipole-expansion theory [15]. They expanded the potentials inside and outside the 
spheres in terms of solid spherical harmonics and gained the values of multipole moments 
Q\m by solving a set of linear equations determined by the boundary conditions at the grain 
surfaces. The upper limit of /, denoted by L, determines how many multipole effects are 
considered: L — 1 refers to the simple dipole approximation and L = oo to the exact cal- 
culation where all multipole effects are included. This multipole expansion theory and our 
DID model are based on different pictures and there exists no direct and exact equivalent 
relationship between them. The calculation of the interparticle energy W,- L using Eq. (6) 
and Clercx et aZ.'s model shows that the DID results reflects, although not exactly, some 
characteristics of the third order(L = 3) multipole-expansion theory, i.e., an octupole effect. 
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This is particularly obvious when 9 = n/2 even in the touching-particle case (See Fig. 5), 
where the point-dipole approximation is known to err considerably. We expect that the 
DID results may be better than the first approximation of Clercx et a/.'s model since the 
former takes higher multipole effects into consideration. Unfortunately the available data 
related with the energy ground state are limited in Ref. [15] and we cannot give a thorough 
comparison. However, from Table 1 in which the ratios of the total induced dipole per par- 
ticle between the bcc and the fee lattice calculated with our model are compared with the 
available results of Clercx et al, we can still see that the DID results is closer to the exact 
solutions than the dipole approximation in a wide range of dielectric mismatch. 



IV. POLYDISPERSE SYSTEMS 

So far we have discussed the interaction energy and ground states in an ER fluid contain- 
ing the same particles. Now we begin to investigate a more "natural" polydisperse system 
in which the particle dielectric constant e p , and consequently the dipolar factor /3, has a 
probability distribution. Let the average interaction energy between a pair of particles in 
this system denoted by W*. Assuming the dipolar factors (3 and (3' of the two particles are 
independent random variables which have the same distribution 

P( a = ^"*-T>- (26) 

we can get W* by taking an average of Eq.(6) over (3 and j3' as 

W* = Wi + W s , (27) 

where Wi is the interaction energy, including the PD results and two DID terms, between 
identical particles with the average dipolar factor (3 , while Wg is the polydispersity correction 

W 5 = -6 f E 2 a 3 {a 2 (3 fi(3cos 2 6 + 1) + u(2(3 2 a 2 + a 4 )(9cos 2 fl - 1)} (28) 

with a and (3$ as the standard deviation and mean of the distribution of (3. On the 
basis of the expression of interparticle energy, we can easily obtain the average self-energy 
of an infinite chain u s * and the interaction energy between two chains u 1 *, in a polydisperse 
system where particle permittivities are distributed randomly, as 

u s *=IF + u s s , (29) 
u t *(p,z)=^(p,z)+u 5 (p,z), (30) 

where the energy increment caused by polydispersity u s s and u s are given by 

u s s = 5 2 uf + (25 2 + 5*)ul (31) 
u\ = 5 2 u T 1 + (25 2 + 5 4 )4 (32) 

In the above equations, 5 is defined as a/(3 to describe the degree of polydispersity, u s , 
u l are the total self-energy and interchain energy for a monodisperse system in which all 
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particles have the same dipolar factor j3 , while uf, u 2 an d u\, u 2 are the corresponding DID 
corrections. They can be calculated using Eq. (10) - (22). 

As a consequence of the change of the chain energies caused by polydispersity, the energy 
gaps between different lattices also vary and hence affect structure formation in ER fluids. 
Fig. 6 shows both Am — (3 and Aug — Po whh 5 — 0.1 to exhibit the effect of polydispersity 
on ground states. Au is defined as u^ t — uj^. when (3q > and u^ t — when (3q < 
(Here "sc" means separated chains) for a system composed of particles (3$, while Aug is that 
caused by polydispersity. It is shown that Aus keeps opposite to Au when f3$ is larger than 
—0.323 and not very near 1.0. In this range the bet lattice has been proved to be the stable 
ground state in a monodisperse system, and the existence of polydispersity tends to destroy 
its stableness. However, the magnitude of the ratio Aug/ Au is only of the order 10~ 2 , so we 
can conclude that, in a polydisperse system where the dielectric mismatch between particles 
and the fluid is not very large, similar structures will form but the ground state is not so 
stable as that in a monodisperse system. 

All the results discussed above are based on the assumption that the dielectric constants 
of particles are distributed at random in the configuration. But the question is: Is the 
distribution really random, especially when the dielectric constants of particles differ greatly 
from each other? In the extreme case, will the phenomenon "phase separation" happen, i.e., 
identical particles gather together upon the application of the electric field? In order to 
investigate this problem, we consider a two-component fluid containing particles with Pi 
and (5 2 . The fractions of particles are supposed to be pi and p 2 , respectively. First assume 
a totally random configuration will form in this system. The average self-energy u s mix and 
interchain energy u l mix in such a system can be calculated by means of the method described 
according to the binary distribution 

P(f3) = Pl 5(f3 - A) + p 2 5((3 - f3 2 ). (33) 

Then consider another configuration in which only identical particles aggregate into "uni- 
form" chains. Two types of chains are formed in this system: one, called chain A, is composed 
of particles and another, chain B, of (5 2 . The average self-energy u s uni f is determined by 

Knif = Pl«* \p=pl +P2U S 1/9=0,, (34) 

and the interaction energy between two uniform chains satisfies 
with 



u. 



l{p, z) = A £ 27^0^) cos(^), (36) 



s=l a 



u. 



„ - P" 1 2n r(f) nF(l 
A is the ratio of /3s of the two chains and other parameters are defined as previously. 
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We have compared the energies in these two configurations and plotted the results in 
Fig. 7. The upper graph in Fig. 7 shows the dependence of Au s = u s mix — u s uni f on f3\ and 
(3 2 , and the lower one shows Au l = u l mix — u l uni j. Note that here u l uni j is the interaction 
between chain A and chain B. The results when we select other chain configuration such as 
chains AA or BB are similar. The fraction of the particles with f3\ is 0.5, and the relative 
coordinates of the interacting chains are chosen as p = \/3a, z = a. It is shown that u s mis 
is always larger than u s uni j, while Au l is usually negative for most values of f3\ and 02 but 
much smaller than Au s in magnitude. Since the self-energy is usually dominant in the total 
interaction energy per particle, chains of identical particles may be more stable than those 
containing different types of particles, particularly when Pi and (3 2 approach the positive 
limit 1.0 and negative limit -0.5 simultaneously. 

The interaction form in polydisperse fluids is more complex than what we have considered 
because it is sensitive to the micro structure formed in the ER solid and the number of 
possible configurations in a polydisperse system is much larger than that in a monodisperse 
one. The energy gap between the random and phase separation configurations may be 
narrowed because of the contribution of interaction between chains. But we can still expect 
that the ground state of such a two-component system may contain quite a few of these 
"uniform" chains instead of a totally random configuration. And the simplicity of our DID 
model will also make it possible to carry out computer simulations in polydisperse systems 
and study the dielectric effects in a more detailed way [23]. 

V. CONCLUSIONS 

In conclusion, we have presented a DID model to study the dielectric effects on structure 
formation in ER fluids. Based on the DID expression of the interacting energy of two 
particles, we have corrected the PD results of the self-energy and the interaction energy of 
chains, including multipole effects partially. Series expressions of these energies are obtained 
and used to calculate the effects of dielectric mismatch between particles and the fluid. Both 
monodisperse and polydisperse cases are discussed and some interesting results are obtained 
as follows. 

(1) In monodisperse systems where particles are negatively polarized, there may exist 
a phase transition form the bet lattice to the configuration of separate chains when f3 < 
-0.323; 

(2) Polydispersity of particle dielectric constants in ER fluids will cause the bet ground 
state not so stable as that in monodisperse systems when the dielectric mismatch between 
the particle and the fluid is not very large; 

(3) When the dielectric constants of particles differ much with each other, identical 
particles tend to aggregate into uniform chains. 

The ER effects in polydisperse system will be studied more deeply in our future simulation 
work based on the DID model. 
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FIGURES 



FIG. 1. The total induced dipole inside a couple of identical spheres A and B upon the appli- 
cation of Eo in the DID model when the dipolar factor (3 > and (3 < 0. 

FIG. 2. (a) two interacting particles; (b) the self-energy u s of an infinite chain; (c) interaction 
energy u l between an infinite chain and a neighboring particle at (p, z). 

FIG. 3. The dependence of u l (p,z) (in units of j^s) on z for different positive (3s in the 
monodisperse system. Solid lines: p = 2a; Dashed lines: p = \/5a; Dotted lines: p = y/Ea. 

2 

FIG. 4. The dependence of u l (p,z) (in units of j^s) on z f° r different negative (3s in the 
monodisperse system. Solid lines: p = 2a; Dashed lines: p = \/5a; Dotted lines: p = \/6a 

FIG. 5. The relative interaction energy (in units of j^s) between a couple of identical particles 
vs the dipolar factor [3 when r = 2 and = tt/2. Solid line: the DID model; Dashed line: the 
multipole expansion theory with L = 3; Dotted line: the multipole expansion theory with L = 1. 

2 

FIG. 6. (a) The energy gap (in units of between the bet and fee lattices Au for a monodis- 
perse system of particles (3o > (solid line) and the corresponding polydispersity correction Au$ 
(dashed lines) ; (b) The energy gap between the bet lattice and separate chains An for a monodis- 
perse system of particles (3o < (solid line) and the corresponding polydispersity correction Aus 
(dashed lines) . 

FIG. 7. The difference of the self-energy Au s (the upper one) and the interchain energy Au 1 
(the lower one) between a random configuration and a system containing only uniform chains. All 
energies are in units of 
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TABLES 



TABLE I. The comparison between the DID model and the full multipole theory proposed 
by Clercx and Bossis . a is defined as e p /e/, where e p and ej are the dielectric constants of the 
particles and the fluid, respectively. A^id is the polarization ratio between the bet and fee lattices 
using the present model. Ai and A^ are from the first-order approximation and the exact results 
of the full multipole theory, respectively. 
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0.984 


0.992 
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a = 10 


1.001 
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